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SUMMARY 

A  single  elastic  fiber  embedded  in  an  infinite  elastic  matrix  is  considered.  An  incident 
plane  SH  wave  is  assumed  in  the  infinite  matrix,  and  an  expression  is  derived  for  the  total 
energy  in  the  fiber  due  to  the  incident  SH  wave.  A  nondimens ional  form  of  the  fiber  energy 
is  plotted  as  a  function  of  the  nondimensional  wavenumber  of  the  SH  wave.  It  is  shown 
that  the  fiber  energy  attains  maximum  values  at  specific  values  of  the  wavenumber  of  the 
incident  wave.  The  results  obtained  here  are  interpreted  in  the  context  of  phenomena 
observed  in  acousto-ultrasonic  experiments  on  fiber  reinforced  composite  materials. 
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INTRODUCTION 


Quantitative  ultrasonic  nondestructive  evaluation  methods  require  analyses  of  wave 
propagation  in  the  material  or  structure  which  is  to  be  evaluated.  In  fiber  reinforced  com¬ 
posite  materials,  wave  propagation  characteristics  depend  on  properties  and  arrangement 
of  the  constituent  materials,  as  well  as  on  the  overall  geometry  and  both  the  microscopic 
and  macroscopic  defect  state  of  the  composite. 

Analyses  of  wave  propagation  in  fiber  reinforced  composites  are  reviewed  in  [l]  and 
[2].  According  to  these  reviews,  wave  propagation  analyses  may  be  divided  into  two  broad 
categories.  In  the  first  category,  the  composite  is  treated  as  an  equivalent  homogeneous 
medium,  and  mechanical  responses  such  as  stress  and  strain  are  described  in  terms  of 
averages  over  the  constitutive  elements.  It  has  been  stated  that  the  homogeneous  or 
equivalent  continuum  approach  is  valid  “when  the  scale  of  the  changes  in  stress  level  (rise 
distance,  wavelength,  etc.)  is  much  larger  than  the  sizes  of  the  constituents  of  the  composite 
(fiber  or  particle  diameter,  fiber  spacing,  ply  spacing,  etc.)”  [l].  In  the  second  category 
of  wave  propagation  analyses,  the  composite  is  modeled  as  a  heterogeneous  mixture  of 
constituents,  and  the  dynamic  interaction  between  the  constituent  materials  is  considered. 
The  heterogeneous  analyses  are  complex,  but  are  necessary  if  the  scale  of  changes  in  stress 
level  is  smaller  than  or  of  the  same  order  as  the  size  of  the  composite  constituents,  or  if 
small-scale  phenomena  such  as  interface  bond  failure  are  of  interest. 

Recent  acousto-ultrasonic  experiments  by  Kautz  [3]  show  some  interesting  wave  prop¬ 
agation  behavior  in  graphite  fiber  reinforced  epoxy  composite  panels.  Kautz  introduces  a 
broadband  signal  into  a  composite  panel,  and  observes  that  the  response  of  the  composite 
appears  to  be  divided  into  a  high  frequency  component  traveling  at  a  relatively  fast  wave 
speed  and  a  low  frequency  component  traveling  at  a  relatively  low  wave  speed.  When  the 
same  broadband  signal  is  introduced  into  a  geometrically  similar  panel  consisting  only  of 
the  epoxy  matrix,  the  response  is  observed  to  consist  only  of  the  low  frequency  component. 
Kautz  hypothesizes  that  the  composite  panel  acts  as  a  heterogeneous  medium,  with  wave 
propagation  along  the  fibers  at  the  higher  wave  speed  and  wave  propagation  through  the 
matrix  at  the  lower  wave  speed. 

In  this  report,  an  infinite  elastic  matrix  containing  a  single  cylindrical  elastic  fiber 
is  considered.  An  incident  plane  SH  wave  is  assumed  to  exist  in  the  elastic  matrix,  and 
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an  expression  is  derived  for  the  total  energy  of  the  fiber  due  to  the  incident  SH  wave. 

It  is  shown  that  the  total  fiber  energy  attains  maximum  values  at  specific  values  of  the 
wavelength  of  the  incident  wave.  The  purpose  of  this  work  is  to  investigate,  in  a  simple 
setting,  the  issue  of  energy  transfer  between  the  matrix  and  the  fibers  in  a  composite 
material.  Recommendations  are  given  for  future  research  which  may  give  further  insight 
into  the  experimentally  observed  phenomena  in  [3]. 

ANALYSIS 

A  cylindrical  fiber  having  circular  cross  section  is  shown  in  Fig.  1.  Rectangular  coordi¬ 
nates  (x,y,z)  and  cylindrical  coordinates  (r,6,z)  are  defined  as  shown  in  Fig.  1.  Displace¬ 
ments  ux,  uy  and  uz  are  defined  in  the  x,  y  and  z  directions,  respectively,  and  displacements 
ur,  ue  and  u2  are  defined  in  the  r,  9  and  z  directions,  respectively.  The  fiber  is  assumed 
to  have  radius  a  and  to  be  of  infinite  extent  in  the  z  direction.  The  fiber  is  assumed  to  be 
elastic,  with  mass  density  pf  and  Lame  constants  /i/  and  A/.  It  is  assumed  that  the  fiber 
is  embedded  in  an  infinite  elastic  matrix  with  mass  density  p  and  Lame  constants  p  and 
A.  Continuity  of  relevant  stresses  and  displacements  is  assumed  at  the  materials  interface 

r  =  a. 

The  infinite  elastic  matrix  is  assumed  to  contain  an  incident  harmonic  plane  wave  of 
the  form 


u^{x,y,z,t)  =  0 
u^{x,y,z,t)  =  0 
u^{x,y,z,t)  =  Re  \ 

[w0e^kx-ut)} 

(1) 

(2) 

(3) 

where  the  superscript  (.')  denotes  the  incident  wave,  u.0  is  the  amplitude  of  the  incident 
wave,  k  is  the  wavenumber  of  the  incident  wave,  w  is  radian  frequency,  t  is  time  and 
j  =  yf-i.  The  wavenumber  k  and  frequency  w  are  related  by 

u  =  ca 

,  k 

(4) 

where 

«-(;) 

1/2 

1  • 

(5) 
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Eqns.  (1),  (2)  and  (3)  define  a  wave  which  propagates  in  the  positive  x  direction  with 
particle  motion  parallel  to  the  fiber  in  the  2  direction.  Such  a  wave,  in  this  context,  is 
denoted  as  an  SH  wave. 

Because  of  the  cylindrical  geometry  of  the  fiber,  it  is  convenient  to  write  the  incident 
SH  wave  in  the  cylindrical  coordinates  (r,  0,  z)  as  [4] 


u^(r,6,z,t)  =  0 

(6) 

u0Hr,d,z,t)  =0 

s 

(7) 

uzHri  zi  t)  —  -Re  <  wo  ^2  enjnJn[kr)  cos  nde  3wt  1 
l  n= 0  J 

(8) 

where  Jn  is  a  Bessel  function  of  the  first  kind  and  of  order  n,  and  where 

fl,  n=0 

1 2,  n=l,2,... 

(9) 

It  is  shown  in  [4]  that  eqns.  (l)  through  (3)  and  eqns.  (6)  through  (8)  are 

mathematically 

identical  representations  of  the  same  incident  wave. 

The  interface  conditions  at  r  —  a  and  the  three-dimensional  equations  of  elasticity 
which  govern  the  infinite  medium  and  the  fiber  are  satisfied  by  the  total  displacement  field 

[4] 

■t 

Jx? 

il 

o 

J 

(10) 

ue(r,0,z,t)  =  0 

►  r  >  a 

(11) 

uz(r,0,z,t)  =  u^\r,0,z,t)  +u [*\r,0,z,t)  > 

(12) 

ur(r,0,z,t)=  0 

(13) 

£ 

<l> 

M 

II 

o 

>  r  <  a 

(14) 

Hoi 

N 

C 

s, 

3 

II 

*40 

<3^ 

tT 

N 

3 

(15) 

where 


u[a\r,0,z,t)  =  Re 
uP(r,0,z,t)  =  Re 


OQ 

Wq  £  AnH^(kr)  cos  nO  e 

n— 0 

°°  1 

—wo  ^2  CnJn{kfr)  cos  nde~Jut  > 

n= 0  J 


(16) 

(17) 
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An  =  (-jn£n/ A)  [/ UfkfJn{ka)J'n{kfa )  -  »kJ'n{ka)Jn{kfa)} 

(18) 

Cn  =  (- jnen/A)fik  [j'n(ka)H^(ka)  -  (fca)Jn(fca)] 

(19) 

A  =  n fk (fca)  J'n[kfa)  -  nkH^' {ka)Jn{kfa) 

(20) 

<*>1  ^ 

II 

(21) 

(nV1* 

c-t=\7f) 

(22) 

The  quantities  An  and  Cn  are  dimensionless  complex  coefficients.  hV  is  a 

of  the  first  kind  and  of  order  »,  and  the  prime  denotes  the  derivative  of 

Hankel  function 

a  function  with 

respect  to  its  argument. 

Eqn.  (16)  represents  a  scattered  wave  which  propagates  outward  from  the  fiber.  As 
r  oo,  the  amplitude  of  the  scattered  wave  approaches  zero,  so  that  far  from  the  fiber 

the  total  displacement  field  approaches  the  displacement  field  of  the  incident  SH  wave. 

Eqn.  (17)  represents  a  refracted  or  transmitted  wave  in  the  fiber.  Eqn.  (17)  is  a 
standing  wave  with  time-harmonic  particle  motion  in  the  z  direction.  Over  any  fiber  cross 
section  *  =  constant,  the  fiber  displacement  in  the  z  direction  varies  as  a  function  of  r  and 

0  according  to  eqn.  (17). 

The  nonzero  fiber  strains  associated  with  the  displacements  given  by  eqns.  (13)  through 

(i5)  are 

1  duin 
lez  ~  r  d9 

(23) 

dvSJ'* 

lrz~  dr 

(24) 

and  the  nonzero  stresses  in  the  fiber  are 

1  dui P 

°0z  -  (*/r  de 

(25) 

3u<'> 

dr 

(26) 

The  elastic  strain  energy  Uf  per  unit  length  of  fiber  in  the  z  direction  is 

/•a  ^ 

Uf=  /  -  {Orzlrz  +  oezlez)  r  dOdr 

Jo  Jo  L 

(27) 
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The  kinetic  energy  Kf  per  unit  length  of  fiber  in  the  2  direction  is 

n2?r  1  .  2 

^f(uif))  rdddr 

The  total  energy  Ef  per  unit  length  of  fiber  is,  therefore, 

Ef  =  Uf  +  Kf 


(28) 


(29) 


The  energy  Ef  given  by  eqn.  (29)  is  a  function  of  time.  In  order  to  obtain  a  single 
numerical  measure  of  energy,  it  is  customary  to  average  the  energy  Ef  over  one  time  period 
of  the  harmonic  oscillation.  The  time-averaged  energy  <  E  >  f  is  defined  by 


where 


1  !T 

>,~tJ0  Efdt 

(30) 

m  27T 

T  =  — 
u 

(31) 

When  eqn.  (17)  and  eqns.  (23)  through  (29)  are  substituted  into  eqn.  (30),  a  long 
algebraic  calculation  gives  the  following  result  for  the  time-averaged  fiber  energy: 

<  E  >,=  |^|g{c*Cni-[(M!{^+1(*/0)  -  Ukfa)Jn+2(hfa) 

+  *^n(^/a)  —  *^n-i(k/a)  Jn+1(k/o)  |  (32) 

\(kfa)  | 


where 


q  _  /  V2.  rc  =  0 

Pn  (1,  n  =  1,2,. .. 


(33) 


and  where  the  asterisk  denotes  the  complex  conjugate.  The  derivation  of  eqn.  (32)  requires 
the  complex  variable  identity 


j(z  +  3‘) . 


(34) 


the  trigonometric  identies  [5] 


1*2% 
Jo  ' 


7 r 


cos  nO  cos  mO  d$  =  —6mn 
0  fin 


(35) 
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sin  nOsinmd 


and  the  Bessel  function  identities  [6] 


={: 


for  n=0  or  m=0 


7 r<5mn  otherwise 


^-Jn{kfr)  =  -kfJn+l(k}r)  +  -Jn{kfr) 
or  ~ 

f  rJ^(kfr)dr  =  —  {J^ik/o)  ~  Jn-i{kfo)Jn+i{kfO,)} 

Jo  L 


j  Jn{kfr)Jn+i{kfr)dr  =  [l  -  J%{kfd)}  -  j  (39) 

fa-Jl{kjr)dr=^\\  +  Jl{ksa)  +  Jl{kfa)-2Y'Jl{k}a)\,  n>  0  (40) 

Jo  r  In  y  k=o  ) 

The  quantity  6mn  in  eqns.  (35)  and  (36)  is  the  Kronecker  delta. 

Eqn.  (32)  is  an  expression  for  the  total  fiber  energy  per  fiber  length  in  the  2  direction. 
In  order  to  nondimensionalize  the  fiber  energy  <E>f,  consider  an  infinite  uniform  elastic 
matrix  containing  only  the  incident  SH  wave  given  by  eqns.  (6)  through  (8),  and  consider 
the  time-averaged  energy  <  E  >m  contained,  per  unit  length  in  the  2  direction,  in  the 
region  0  <  r  <  a.  By  a  calculation  similar  to  the  calculation  leading  to  eqn.  (32),  the 
time-averaged  energy  <  E  >m  is 


<  E  >m=  o 

4 


jr^D*nDnj-  ^j-{jn+i{ko)  ~  Jn{ka)Jn+2{ka) 

+J*{ka)  -  Jn_i{ka)Jn+i(ka )  j  +  nJ^(ka) 


where 


Dn  —  triJ 


A  nondimens  ional  fiber  energy  <  e  >f  is  now  defined  by 


<  e  >f= 


<E_>2 

<  E  >m 


The  dimensionless  quantity  <  e  >/  gives  a  measure  of  the  fiber  energy  relative  to  the 
energy  the  matrix  would  contain  if  there  were  no  fiber  present. 
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By  introducing  the  nondimensional  parameters 


eqns.  (43)  and  (19)  can  be  written  as 

<  E  >f 

v  - —  J 


<e  >f= 


<  E  >, 


y°°  C*C  —  /  Jn+A'1ka)  ~  Jn(lka)Jn+2{lka)  1  .  np(~un\ 

=  r  ^n=°  "  [  2  l  +Jfaka)  -  J  + 

Y°°  D*D  —  [(fc°)g  /  Jn+l{ka)  -  Jn{ka)Jn+2{kci)  \  I  n  Tiflcn) 
2.n=0  2  j  +  j2(fca)  -  Jn_,  (ka)  Jn+1  (ka)  J  +  ^ ^ 

i-jnen)  [j^(ka)H{nl)(ka)  -  H{nlY (ka)Jn(ka)'\ 


rfl')Hi1\ka)J^ka)  -  H^iy (ka)  Jn(lka) 


-  rP 


respectively,  where 


From  eqns.  (46),  (47)  and  (48)  it  can  be  seen  that  the  nondimensional  fiber  energy  <  e  >/ 
depends  only  on  the  nondimensional  parameters  ka,  rp  and  rp.  The  parameter  ka  is  a 
nondimensional  wavenumber  of  the  incident  wave,  and  can  be  written  in  terms  of  the 
incident  wavelength  A,-  as 

.  2tt  a 

ka  —  — — 

At 

or  in  terms  of  frequency  as 


The  nondimensional  fiber  energy  <  e  >/  is  plotted  as  a  function  of  the  nondimensional 
wavenumber  ka  in  Fig.  2  for  the  numerical  values  rp  =  2.0  and  rp  =  25.0.  These  numerical 
values  correspond  to  an  E-glass  fiber  in  a  PMR  (polyimide  monomeric  reactant)  matrix 
[7].  In  Fig.  3,  the  nondimensional  fiber  energy  <  e  >/  is  plotted  as  a  function  of  ka  for  the 
numerical  values  rp  =  1.5  and  rp  =  7.5.  These  numerical  values  correspond  to  a  Celion 
6000  fiber  in  a  PMR  matrix  [7j.  The  numerical  values  used  to  generate  Fig.  3  are  estimates 
of  the  material  properties  of  the  graphite/epoxy  composite  used  in  [3].  The  infinite  series 
in  eqn.  (46)  are  computed  by  truncating  when  the  absolute  value  of  the  individual  terms 
becomes  less  than  10-12  times  the  absolute  value  of  the  accumulated  sum. 


DISCUSSION 


In  both  Fig.  2  and  Fig.  3,  the  nondimensional  fiber  energy  <  e  >/  is  approximately 
equal  to  one  when  ka  is  close  to  zero,  or,  equivalently,  when  the  wavelength  of  the  incident 
wave  is  long  compared  to  the  fiber  diameter.  This  means  that  for  long  wavelengths,  the 
energy  in  the  fiber  is  approximately  the  same  as  the  energy  which  would  exist  in  a  uniform 
matrix.  In  other  words,  the  presence  of  the  fiber  does  not  significantly  affect  the  energy 
distribution  for  long  wavelengths.  As  ka  becomes  larger  (or,  equivalently,  as  the  wavelength 
becomes  smaller)  the  energy  in  the  scattered  wave  becomes  larger,  and  the  energy  in  the 
fiber  becomes  smaller  than  the  energy  which  would  exist  in  the  region  0  <  r  <  a  in  the 
uniform  matrix. 

The  most  interesting  features  of  Fig.  2  and  Fig.  3  are  the  local  maxima  in  the  values 
of  <  e  >/.  For  an  incident  wave  containing  a  broad  band  of  wavenumber,  the  fiber  energy 
will  be  greatest  at  frequencies  corresponding  to  the  values  of  ka  where  the  local  maxima 
occur.  If  the  fiber  response  could  be  measured,  the  local  maxima  would  appear  in  a 
frequency  spectrum  of  the  measured  response.  Thus  it  is  sensible,  at  least  in  this  simple 
case,  to  suppose  that  certain  frequency  components  in  a  measured  response  may  contain 
information  about  the  fiber. 

The  incident  SH  wave  assumed  here  is  independent  of  the  spatial  coordinate  z  and  so, 
therefore,  is  the  displacement  field  in  the  fiber.  Therefore,  it  is  impossible,  in  this  analysis, 
to  discuss  the  issue  of  wave  propagation  along  the  fiber  in  the  z  direction.  The  analysis  of 
an  incident  wave  which  has  a  z  component  may  allow  consideration  of  wave  propagation 
along  the  fiber.  Equations  which  describe  the  scattering  and  refraction  of  such  an  obliquely 
incident  wave  are  given  in  [4]. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  report  is  to  investigate  the  transfer  of  energy  to  a  single  fiber  in 
an  infinite  elastic  matrix  containing  an  incident  SH  wave.  An  expression  is  derived  for  the 
energy  in  the  fiber,  and  it  is  shown  that  the  fiber  energy  attains  local  maxima  at  specific 
values  of  the  wavenumber  of  the  incident  wave. 

Future  research  should  consider  the  problems  of  incident  P  and  SV  waves.  As  men¬ 
tioned  above,  an  obliquely  incident  wave  may  allow  analysis  of  wave  propagation  along 
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the  fiber.  Also,  properties  of  the  scattered  wave  given  by  eqn.  (16)  should  be  investigated, 
since  the  scattered  wave  may  be  available  for  experimental  measurement.  It  would  be  in¬ 
teresting  to  see  if  the  energy  flux  of  the  scattered  wave  has  wavelength  dependence  similar 
to  the  wavelength  dependence  of  the  energy  of  the  fiber. 
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Fig.  1  Elastic  fiber  embedded  in  infinite  elastic  matrix  containing  incident  SH  wave 
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